Dark Energy Models in the w — w' Plane 
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We examine the behavior of dark energy models in the plane defined by w (the equation of state 
parameter for the dark energy) and w' (the derivative of w with respect to the logarithm of the 
scale factor). For non-phantom barotropic fluids with positive squared sound speed, we find that 
w' < 3w(w + 1), the opposite of the bound on quintessence models previously derived by Caldwell 
and Linder. Thus, these barotropic models and quintessence models for the dark energy occupy 
disjoint regions in the w — w' plane. We also derive two new bounds for quintessence models in the 
w — w' plane: the first is a general bound for any scalar field with a monotonic potential, while the 
second improves on the Caldwell-Linder bound for tracker quintessence models. Observationally 
distinguishing barotropic models from quintessence models requires cr(w') < 1 -+- w. 



I. INTRODUCTION 

Observations indicate that roughly 70% of the total en- 
ergy density in the universe is in the form of dark energy, 
which has negative pressure and drives the late-time ac- 
celeration of the universe (seeRef. Q for a recent review, 
and references therein). 

The simplest model for this dark energy is a cosmo- 
logical constant, but many other alternatives have been 
put forward. One possibility, dubbed quintessence, is a 
model in which the dark energy arises from a scalar field 
El [J 0- E| In another class of models, the dark en- 
ergy is simply taken to be barotropic fluid, in which the 
pressure p and energy density p are related by 



P = /(p). 



(1) 



The prototype for this sort of model is the Chaplygin 
gas 0, 13 , which can be considered a special case of the 
generalized Chaplygin gas |9(. Although originally pro- 
posed as unified models for dark matter and dark energy, 
the Chaplygin gas and generalized Chaplygin gas mod- 
els have also been examined as models for dark energy 
alone [HI [HI [2, [l|, H3 ■ Interest in barotropic fluids 
as dark energy has not been confined to the generalized 
Chaplygin gas. Many other models have been proposed, 
including, e.g., the Van der Waals model 0] and the wet 
dark fluid model |l6j. 

Obviously, it is important to determine whether these 
models can be distinguished by their observational conse- 
quences. In this vein, Caldwell and Linder |17| examined 
the behavior of quintessence models in the plane defined 
by the quantities w and w' for the dark energy. Here w 
is the ratio of pressure to density for the dark energy: 



w 



p/p, 



(2) 



and w' is the derivative of w with respect to the logarithm 
of the scale factor a: 



w 



dw 
dln(a) 



(3) 



Caldwell and Linder showed that quintessence models 
in which the scalar field potential asymptotically ap- 



proaches zero can be divided naturally into two cat- 
egories, which they dubbed "freezing" and "thawing" 
models, with quite different behavior in the w — w' plane. 

Here we extend and generalize these results. In the 
next section, we determine the behavior of general non- 
phantom barotropic fluids with positive squared sound 
speed in the w — w' plane, and we show that these mod- 
els occupy a region of phase space which is disjoint from 
that occupied by quintessence models. In Section 3, we 
re-examine the lower bound on w' from Ref. |17| ; we de- 
rive a new bound for a more general class of quintessence 
models, and sharpen the bound from Ref. [l7j for the 
case of tracker potentials. These results are discussed in 
Section 4. Our results support the argument that the be- 
havior of dark energy models in the w — w' plane provides 
a useful discriminator for such models. 



II. BAROTROPIC FLUIDS 

As noted in the previous section, a barotropic fluid is 
one for which the pressure is purely a function of the 
density, p — f(p). For such models, we have 



, dw dp 

w = 



dp din (a) 

The factors in this equation are given by 



and 



dw 
dp 

dp 



1 /dp 
p \dp 



-3(l + w)p. 



dln(a) 

Then the expression for w' becomes 



w 



(4) 



(5) 



(6) 



(7) 



In order to derive useful constraints from equation J7J) , we 
need to assume something about the behavior of dp/ dp 
and w. We will assume first that we are not dealing with 
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phantom models, so that 1 + w > 0. Second, we note 
that dp/dp = Cg, and we will assume further that c 2 > 0, 
to avoid instabilities in perturbation growth. These con- 
straints are satisfied by most barotropic fluids of interest, 
although it is also possible to construct models which vi- 
olate both assumptions (see, e.g., Ref. 0]). With these 
assumptions, equation (JJJ gives the limit 



w' < 3w(l + w). 



(8) 



This limit is shown in Fig. 1, where have used the same 
scale as Ref. for easy comparison. Note that this is 
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FIG. 1: Solid curve gives the upper bound on w' as a function 
of w for non-phantom barotropic models with positive squared 
sound speed. Dotted curves give w' as a function of w for 
Chaplygin gas dark energy with (top to bottom) a = 0.5, 
a = 1, a = 2, while a = lies on the solid curve. Dashed 
curves give w' as a function of w for the wet dark fluid model 
with (top to bottom) to, = 0.1, w* = 0.2, tu, — 0.3. 

the exact opposite of the limit derived in Ref. 0] for 
"freezing" quintessence models: 

w' > 3w(l + w). (9) 

Thus, barotropic fluids and freezing quintessence models 
occupy non-overlapping regions in the w — w' plane. The 
bound in equation ||SJ| is saturated when c 2 = dp /dp = 0, 
which corresponds to the case p = constant. In this case, 
we have w' = 3u>(l + w). Models of this kind are partic- 
ularly interesting, since they correspond to a barotropic 
fluid that behaves exactly like a mixture of dark mat- 
ter (dust) and a cosmological constant. Such models can 
be constructed in the context of either the generalized 
Chaplygin gas 18] or fc-essence |T^ . 

Now consider several illustrative special cases. For the 
generalized Chaplygin gas, the equation of state is given 
by 

p=-4 ( io ) 



where A and a are constants. (The Chaplygin gas corre- 
sponds to the special case a = 1). Then it is straightfor- 
ward to use equation J2J to derive 



to' = 3(1 + a)w(l + w). 



(11) 



This behavior is displayed in Fig. 1 for several values of 
a. The a — case saturates the upper bound on w' , 
as it corresponds to the previously-noted p = constant 
model, while the original Chaplygin gas model (a = 1) 
gives w' = 6w(l + w). 

Another barotropic model is the wet dark fluid model 
of Ref. [lj| , with equation of state 



p = w*(p - p„), 



(12) 



where u>* and are constants. Again, it is straightfor- 
ward to use equation J7J to derive the result 



w' = —3(1 + w)(w* — w), 



(13) 



which is displayed in Fig. 1 for several values of w*. 

While we have assumed that the fluid characterized by 
w' and w acts strictly as dark energy, all of the results 
derived here can also be applied to unified dark energy 
models, in which the barotropic fluid serves as both dark 
matter and dark energy (indeed, this is one of the main 
motivations for the Chaplygin gas model and its vari- 
ants). However, in the case of unified dark energy, w and 
w' refer to the entire fluid, not just to the dark energy 
component. 



III. QUINTESSENCE 



= 0, 



The equation of motion for the quintessence field 

where the Hubble parameter H is 

a /px 1 / 2 
a ~ V3/ 



H 



is 



(14) 



(15) 



and p is the total density. (We assume an £1 = 1 universe 
and take 8nG = 1 throughout). Since we are interested 
in the late-time evolution of the dark energy, we assume 
that p = pm + p^, where pu is the density of the matter, 
scaling as pm oc a -3 , and p^ is the scalar field energy 
density, given by 



P<l> 



The pressure of the scalar field is 



(16) 



(17) 



Caldwell and Linder [l7j examined models in which 4> 
evolves toward the state V(<fr) = 0. They distinguished 
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two cases: "thawing" models, in which w ~ — 1 initially, 
but w increases as <fi rolls down the potential, and "freez- 
ing" models, in which w > — 1 initially, but w approaches 
— 1 as the field rolls down the potential. The bounds de- 
rived in Ref. |17| are based on a combination of plau- 
sibility arguments and numerical simulations. For the 
"thawing" models, these limits are 

1 + w < w' < 3(1 + to), (18) 

while for the "freezing" models, 

3to(l + to) < to' < 0.2w(l + w). (19) 

As noted in the previous section, the allowed region for 
barotropic fluids in the w — w' plane lies adjacent to the 
lower bound on w' for quintessence models proposed in 
Ref. yjj- Therefore, it is the lower bound in equation 
(|19fl which is of greatest interest here. 

We first generalize the Caldwell-Linder lower bound 
to a wider range of potentials. Consider a scalar field 
evolving in an arbitrary monotonic potential V(4>), not 
necessarily asymptotically zero. We will take V'(4>) < 0, 
but of course the argument is identical in the opposite 
case. Now, following Steinhardt et al. @, we note that 
the equation of motion for (f> can be rewritten as 



3(1 + w) 

v ~ V a 



<t> 



X' 

1 + T 



where X is defined as 

X = In 



1 



to 



(20) 



(21) 



and X' is the derivative of X with respect to ln(a), so 
that X' and w' are related via 



X' 



2w' 



(1 - to)(l + to) ' 



(22) 



Now we note that for a scalar field rolling downhill in 
a monotonically decreasing potential, the left-hand side 
of equation H20|) is always positive. This, in turn, implies 
that 1 + X' /6 > 0, which translates into the following 
bound on to': 



> -3(1 - tu)(l +w). 



(23) 



While not as tight as the Caldwell-Linder bound, equa- 
tion l|23|) applies to a more general class of models; it as- 
sumes nothing about tracking behavior or "freezing" of 
the scalar field, nor is it confined to potentials for which 
V(4>) — > asymptotically. This bound is saturated by 
a scalar field rolling on the constant potential V = Vq. 
For this case, uu evolves from to~+ltoto« — 1, with 
to' = -3(1 - to) (1 + to) [13. The bound inequation [|23)l. 
along with the lower bound on to' from Ref. ^t|> are 
shown in Fig. 2. 

We can tighten this bound considerably for the case 
of "tracker" models, i.e., models for which w is initially 
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FIG. 2: Solid curve gives the lower bound on w' as a function 
of w for "freezing" quintessence models from Ref. fl^l . Dotted 
curve gives our lower bound on w' as a function of w for 
quintessence models with an arbitrary monotonic potential 
(equation 1231 . Dashed curve gives our lower bound on w' as 
a function of w for tracker quintessence models (equation l3H . 



constant when <C pm, but which evolve toward w = 
— 1 when w pm- Again, following Ref. 0, we define 
r to be given by 



r = 



V"V 



yr. 



(24) 



If u>b is the equation of state parameter for the back- 
ground fluid, then tracking behavior with w < wb occurs 
when r is roughly constant with r > 1 ||. Note that T 
is exactly constant for power law and exponential poten- 
tials. In particular, for the potential V oc <p~ a , we have 
r = 1 + 1/a. Taking the derivative of equation (|20l) with 
respect to <j> gives 



2X" 



1 



X' 



(1 + to) (6 + X') 2 
3(wb — w) 1 — n<£ 



2(1 + w) 6 + X' 



1 



6 + X' 



(25) 



Note that this equation differs from the corresponding 
equation in Ref. [(| because we do not take the limit 
— ► 0; all of the interesting evolution occurs when p^ 
becomes significant compared to pm- It does however, 
agree with the expression previously derived in Ref. |2lj | , 
and in the limit where — > 0, equation (1251) reduces to 
the corresponding equation in Ref. [3 . 

When the tracker solution is reached, and the back- 
ground fluid dominates, w is a constant. We will be 
interested in the late-time evolution of the quintessence 
field, so we assume that the background fluid is cold dark 
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matter, with wb = 0. Then for the tracker solution, tak- 
ing X" = X' = and f2^ = in equation (|25[) gives 



-2(r-i) 
2(r-i) + i 



(26) 



In this regime, increases with time; when fi^ is no 
longer negligible with respect to 1, the tracker solution no 
longer applies. At this point, w decreases, asymptotically 
approaching w = — 1. (See Ref. 0| for a discussion of 
the evolution in this regime). 



0. Initially, X' = 0, but as fi^ becomes nonnegligible, 
X' decreases. It reaches some minimum value, X' min , 
but then increases back to a value near zero as Q<f, — > 

1. This behavior is illustrated in Fig. 3, in which we 
show the evolution for X' as a function of Vt^ for several 
power-law potentials. The value of X' min is the key to our 



on X' min . We first note that < 1; applying this limit 
(with w < 0) to equation i|27[l yields 



X' 



2(r-i)(i 



w 



"" n> 6 (l-™)+2(r -!)(! + «;)• 



(28) 



Since w' < for these tracker solutions, we have X' < on X' 



The right-hand side of equation l|28l) increases with de- 
creasing w, so a lower bound on the right-hand side is 
achieved by assigning the maximum allowed value to w. 
In tracker models, w decreases from the initial value given 
by equation (|26|l . so substituting the value for w given 
by equation 1)26(1 into equation l(28|l gives a lower bound 



i2(r-i) 
mm 6(r -!) + !' 



(29) 



Taking X' > X' min , and using equation l(2^|) to convert 
from X' to w' we obtain the limit 
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FIG. 3: The evolution of X' as a function of Q$ for the po- 
tential V(4>) oc 4>~ a with a — 2 (dashed curve), a = 1 (solid 
curve) and a — 0.2 (dotted curve). 



w > 



6(r-i) 
6(r-i) + i 



(1 - w)(l + w). 



(30) 



This is obviously a model-dependent result, since it de- 
pends on the value of T for which tracking behavior 
is achieved. However, we can use it to derive a more 
general limit by noting that, since T > 1, we have 
6(r - 1)/[6(T - 1) + 1] < 1, giving 



w' > -(1 - w)(l + w). 



(31) 



Equation (|3 II) gives a lower bound on w' for tracker mod- 
els. This bound is shown in Fig. 2. This is clearly a very 
conservative limit; for example, it is obvious from Fig. 3 
that the value of X' min is often obtained for a value of 
fi^ larger than is consistent with observations. However, 
equation (I31|l does give a tighter bound than that derived 
in Ref. , further separating quintessence models from 
the barotropic fluid models discussed in the previous sec- 
tion. 



IV. DISCUSSION 



calculation, since it gives a lower bound on w' through 
equation 1221) . 

To find this minimum value for X 1 , we simply take 
X" = in equation Q25p. where we assume also that 
wb = and that T is nearly constant. This gives 

x , = (i-^V + 2(r-i)(i + w) 

{l - w) + 2{T -!){! + w) ' 1 ' 

On the right-hand side of this equation, w, fi^, and T 
are to be evaluated at the point at which X' — X' min . 
In general Qa, and w evolve in a complex fashion (see, 
e.g., Ref. 122]), so equation lf77|l does not have a simple 
solution. However, since we seek a lower bound on X', 
it is sufficient to use equation l|27|l to find a lower bound 



Our results illustrate the usefulness of the w — w' plane 
as a means for observationally distinguishing dark en- 
ergy models. Non-phantom barotropic fluid models with 
positive squared sound speed and "tracker" quintessence 
models occupy disjoint regions in this plane. Even our 
less stringent general limit on quintessence models with 
monotonic potentials is inconsistent with a wide class of 
barotropic fluid models, as can be seen by comparing 
Figs. 1 and 2. 

Supernova la measurements have narrowed the 

range of possible values for w; in particular, they are 
consistent with w » —1. The value of w' is less well 
constrained (see, e.g., Ref. [H). The results of Ref. [13 
indicate that the resolution on a measurement of w' must 
be of order 1 + w in order to distinguish "freezing" and 
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"thawing" quintessence models. We reach a similar con- 
clusion with regard to barotropic fluids and quintessence 
models: observationally distinguishing between these two 
classes of models also requires that a(w') < 1 + w. 

Finally, we note that there is a relation between the 
w — w' parametrization of Ref. 0] and the "statefinder" 
pair introduced in Ref. [2{| . The statefinders r and s are 
defined by 



a 



r - 1 
3(9-1/2)' 



(32) 
(33) 



where q = —ad /a 2 is the deceleration parameter. Then 
r and s are related to w and w' via 25] 

9 3 
r = 1 + -Q, DE w{l + w) - -n DE w', (34) 



s = 1 + w- — , 
3w 



(35) 



where Qde is the fraction of the total energy density 
contributed by the dark energy. Note, however, that the 
mapping between u>, w' and r, s is non-trivial because of 
the factor of Qde hi equation (|34|) . For non-constant w, 
the evolution of £Ide depends in a complicated way on 
both a and w(a). Hence, the statefinder parameters en- 
code somewhat different information about the behavior 
of the dark energy than do w and w' . The statefinder 
parameters are more directly related to observable quan- 
tities, while the w — w' parametrization describes more 
directly the physical properties of the dark energy. 
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